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I. INTRODUCTION 

The effects of parity violation in tlie two-level approximation for the quantum dynamics of a pair of strictly isolated 
molecular enantiomers [ij was considered some years ago by Harris and Stodolsky (HS) [Ij. Based on this model, 
they pointed out the interesting possibility that tunneling should exist leading to oscillations in the optical activity 
(OA). They argued moreover that such a system should be sensitive to extremely small energies and could be used 
to observe the presence of the weak interaction via oscillations of the optical activity about a non-zero value. Their 
basic idea continues to motivate detailed proposals for measuring the optical activity a highly nontrivial pursuit. 
In this vein, it is also worth mentioning the many varied independent experimental efforts and proposals aimed at 
detecting parityviolation (PV) in chiral molecules. These techniques include vibrational- rotational Q, electronic [1- 
0], Mossbauer Q, and nuclear magnetic resonance (NMR) spectroscopy Q. Proposals for measuring the PV energy 
difference in crystallization [10] and in solubility experiments [IjiJ have also been considered. To date, however, no 
effects of parity violation in chiral molecules have been experimentally observed. The challenge has been taken up 
recently by a multi-disciplinary consortium to employ high-resolution laser spectroscopy for a first observation of PV 
in chiral molecules [l^ . 

Returning to the HS scheme, they also recognized that obstacles to the observation and detection of these chiral 
oscillations in the OA will come from interactions with the surrounding medium such as collision effects, thus tending 
to induce relaxation phenomena [l3j . Radiative processes, fundamental interactions with the radiation field, are also 
important and their inclusion necessarily implies adopting a complex multi-level treatment [l4| . Indeed, the effects 
of both coUisional and radiative processes can be simulated approximately by phenomenologically adding complex 
energies Ek — fPfc, to the hamiltonian spectrum, where Tk denotes a decay width [l^. Apart from effects of collisions 
between the enantiomers and a background medium (i.e., gas, liquid or solid), the double well model implies a tower of 
excited electronic- vibrational states for the chiral molecules [16,], and transitions between these states, induced by an 
external field, is expected to modify the fundamental oscillation period and amplitude of the OA of the ground-state 
enantiomers. These considerations would seem to invalidate the simplest two-state HS description. It is clear that 
the multi-state nature of real molecules should be taken into account in order to obtain a more accurate description 

Motivated by the above-mentioned considerations, we probe somewhat further in the quantum-mechanical model 
of HS. Therefore our objective in this paper is to consider a multi-state approach applied to the racemization problem 
and demonstrate that an effective two-level Hamiltonian description can be derived from perturbation theory, thus 
legitimizing this description. We derive the explicit form of the energy spectrum of the reduced two-level system 
in terms of the appropriate interaction hamiltonian. Our aim is to keep the discussion as general as possible, thus 
our basic assumption is the existence of a hermitian interaction hamiltonian responsible for inducing the transitions 
from the chiral molecules to higher (e.g., electronic, vibrational, and/or rotational) levels. Given this interaction, 
we then adapt straightforwardly the Weisskopf-Wigner (WW) approximation scheme [13] , originally used for solving 
the line width problem in atomic transitions jisj . to the case of chiral molecules treated as a complex multi-state 
quantum mechanical system. The effect of transitions is two-fold: on the one hand, they lead to corrections to 
the enantiomer mass matrix, thus lifting the enantiomeric mass degeneracy and hence yielding corrections to the 
oscillation period; and on the other, they also in principle allow for effects of decay channels of the enantiomers. 
However, in the context of real molecules, such energy-conserving decay channels are most likely absent, as these 
would entail the fragmentation or dissociation of the enantiomers into smaller molecular units. This notwithstanding, 
the WW approach convincingly establishes a formal analogy between the dynamics of the racemization of molecular 
enantiomers and the interference and decay effects predicted [l^ and observed [1^ [2l| in the neutral kaon system of 
elementary particle physics. These are both examples of multiple-level quantum systems with oscillations, although 
the underlying physics in each system (chiral molecules, JT-mesons) is of course radically different. 

We give a brief overview of the simplest Harris-Stodolsky two-state approach in Sec. [Hi Motivated by the multi-level 
nature of real molecules pointed out above, we approach this problem by applying the Weisskopf-Wigner approximation 
scheme in Sec. IIIII This permits one to derive an effective two-level description which takes into account transitions to 
an arbitrary number of multiple excited states. The "end" result is expressed as a compact 2x2 mass matrix which 
acts in the space of the two enantiomers. The associated eigenvalues and eigenvectors depend on whether one assumes 
the putative interaction hamiltonian is CPT or T invariant, and these invariance principles also effect the final form 
of the OA and the oscillation period. Our goal here is to derive an effective Hamiltonian description for the transitions 
to and from multiple levels, not a study per se of the effects that PV may itself lead to (2^ - l29| . although these can be 
included in princi ple as part of the overall perturbation. Dynamical analogies between enantiomers and the neutral 
ii'-meson system [ig, |30, [3l[ , are discussed in Sec IIVI We comment on the novelty of the approach developed in 
the Discussion, and contrast it briefly to other current theoretical works. The WW method is summarized in the 
Appendix. 
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II. HARRIS-STODOLSKY TWO-STATE MODEL 

The Hamiltonian H describing the dynamics of a pair of strictly isolated chiral molecules |L >, |i? > is f3|: 

H = Eol + 5a^ + ea, ^ Ho + ea^, (1) 

where S =< L\Ho\R > is a parameter related to barrier height and e is the energy shift due to parity violation. 
For e = 0, the eigenvectors |+ >= {\L > +\R >)j\/2 and |— >= (|L > —\R >)/^ are states of definite parity: 
P\+ >= 1+ > and P\— >= — |— >. The eigenvalues are given by E± = Eo±S, and level splitting is due to tunneling 
alone. The system wave function obeys — |^ >, and in terms of the parity basis, the general time-dependent 

solution is given by 

ae-*(^«+'^)*/'*|+ > +6e-'(^«-'^)*/^|- > (2) 

where the values of a, b incorporate the initial condition |^'(0) >. Thus, if we prepare the system to be initially in the 
chiral state \L > , then a = b = l/\/2, and at any later time t >0 the wavefunction is given by 

\-qj{t) >= e-'^"*/''{cos{6t/h)\L > -ism{St/h)\R > ). (3) 

Then PL(t) = | < > P = cos^i^) and PR{t) = | < > P = sin^(f ) are the probabilities for the system 

to be in chiral state L or to make a transition to the state R, at time t, respectively. Note that -Pl(0 + Pnit) — 1- 
The optical activity is given by Q 

e(t) = e™ax(PL(i)-Pfl(i)) 

,2St. 

= Qmax COS [—). (4) 

For parity violation, |e| > 0, and the eigenvalues are i?i_2 ~ Eq± {S"^ +e^)^/^ and the associated eigenstates are given 
by l^'i >= cos(j)\L > +sin0|i? > and 1^*2 >= — sin0|L > +cos(f)\R >, with mixing angle defined via cot 20 = | [l]. 
If the system is initially prepared to be in the chiral state \L >, then at any later time its wavefunction is given by 

>= e-*^°*/'*(^(cos2(/>e-*'^*/'"' + sin2 0e+'^*/'*)|L > +cos(^sin0(e-'^*/'* - e'^*/'*)|i? > ), (5) 
where A = (S^ + e^Y^"^. The probability to make a transition to a state \R >, is therefore given by 



PR{t) = I < > 



|2 



= 4 sin^ (/) cos^ <j) sin^ [\/ + dH/fi) 
52 



e2 + ,52 

The probability to be in the state \L > is: 



in^ {^/7^^+^t/n). (6) 



Pdt)^\<L\^{t)> 



= cos2( — ) + cos"' 20 sm2( — ), 
Ti n 

= cos2(_) + ^sm2(_). (7) 
Once again, we have Phit) + PB.{t) = 1. The corresponding optical activity (OA) is then calculated to be 

This latter formula is the starting point for recent proposals for measuring the parity violating energy difference e 
between enantiomers Q. 



4 



III. MULTI-LEVEL PROBLEM: TRANSITIONS TO MANY STATES 



We start with mass degenerate pair \L >, \R > of chiral enantiomers: vriL — = to, considered as ground states 
of a Hamiltonian H^. We then include transitions from these states to other eigenstates or levels of Hq induced by 
external fields (radiation, thermal effects, etc.). The problem we thus consider is the time evolution of a state initially 
prepared as a superposition of |L > and \R >. We regard the potential well barrier, the parity- violating energy 
difference e as part of the overall perturbation Hi. So we decompose Hi = Sa^ + ecr^ + H2, where the first two terms 
act only within the two-dimensional subspace of the ground-state enantiomers, and H2 induces transitions from these 
to the other (electronic-vibrational) levels. These specific considerations show up only at the stage where we display 
the explicit matrix elements of the effective two-level hamiltonian. The following analysis is, however, independent of 
the explicit form of the overall perturbation Hi. 



A. Mass matrix 



We have a system described by the Schrodinger wave function |^'(t) > whose time evolution is given by (h = 1) 

«^|*(<) >=(^o + ^i)|*(<) > . (9) 

We write the full Hamiltonian H = Hq + Hi, here \L > and |_R > are the two degenerate discrete eigenstates of Hq, 
that is, a pair of mirror-image enantiomers, and the perturbation Hi induces transitions from these states to other 
(possibly unbound) eigenstates |fc > of Hq and possibly also between \L > and \R >. The problem to solve is the 
time evolution of a state initially prepared as a superposition of the two degenerate chiral states. In terms of the 
interaction representation, the state-vector 

\ij{t)>=e'"«'\^{t)> (10) 

satisfies the equation 

i4|V(t) >= e''^°*i?ie-''^°*|V(i) >= H'i\ij{t) >, (11) 
at 

so that the time dependence of \ipit) > arises solely from the perturbation term Hi. We expand the interaction 
representation wave function in terms of the complete set of eigenstates: 

\^it) >^a{t)\L > +b{t)\R> +Y,Ckit)\k >, (12) 

k 

subject to the initial conditions a(0) = ao,6(0) = bo, and Cfc(O) = 0. Then Eq. (|1H) leads to the following set of 
coupled equations for the probability amplitudes: 

.da^ ^ ^ ^^^^^ ^ ^(^^^ ^ L|i/i|i? > b{t) + L\H[\k > Ck{t), (13) 

k 



.db{t) 
dt 



= < R\Hi\L > a{t)+ < R\Hi\R > bit) +J2< R\H'i\k > Cfc(i), (14) 



= < k\H'i\L > a{t)+ < k\H'i\R > b{t) + ^ < k\H'i\j > c,{t). (15) 

j 

We may omit the prime on Hi for the matrix elements taken within the \L >, \R > subspace in Eqs. and dTi]) , 
since these are assumed to be degenerate in mass. Applying the Weisskopf-Wigner [17] approximation procedure to 
this multi-level system yields the effective two-level quantum-mechanical description (see Eq. IA14I) 

i^* = fM-ir)*. (16) 



dt 

The hermitian mass and decay matrices M ~ and F = have the explicit matrix elements at order 0{Hi) 

^ Ek-m 



5 



and 

Tap = 2tt^< a\Hi\k >< fc|Fi|/3 > S{Ek - m), (18) 

k 

where the indices a, (3 stand for the states L or R. 

The only assumptions that go into obtaining the resuh in Eqs. (|16p - (|18p are that the dynamics is determined by 
the time-dependent Schrodinger equation, the higher order terms 0{Hf) are neglected, and that the Hamiltonian H 
is hermitian (see the Appendix for details). This result allows for the possibility of energy-conserving decay channels 
through the decay matrix T. Assuming that the chiral enantiomers are not unstable, there will be no decay, that 
is, no fragmentation nor dissociation of the enantiomers into other molecular species. Barring this possibility, there 
will be no contribution from Eq. (fT8|) because Ek > m. One might nevertheless be tempted to think that T could 
automatically account for effects of elastic collisions. Indeed, the full operator structure M— zF of the right hand side of 
Eq. p6p implies complex energies for the spectrum, and is reminiscent of the terms that are added phenomenologically 
to the molecular hamiltonian as a way of simulating approximately the effects of collisions and radiative effects [ij] . 
However, its inclusion would imply exponential decay in the probabilities PL(t), Pji{t) themselves, as well of course 
in the optical activity ^ (Pr — Pr), so that e.g., PL{t), Pnit) ^ as time increases, whereas coUisional effects yield 
instead Pl(^), Pnit) ^ ^ I2j]- Henceforth, we set F = in the remainder of this paper. 



B. Eigenvalues and eigenvectors: CPT and T-invariance 

The underlying assumed invariance affects the form of the eigenvalues and eigenvectors of Eq. (jl6p and it is therefore 
of interest to consider independently the implications of first CPT and then T invariance. The former applies only 
for the case of the so-called CP-enantiomers, which require the existence of mirror image molecules composed of 
antiparticles that is, the CP partner of |L > is \R >, the anti-right handed molecule, whereas the CP partner of 
|i? > is |L >, the anti-left handed molecule. 

Thus the solution of the eigenvalue problem 

M*±==A±*±, (19) 
assuming CPT invariance, so that A/n — M22 (and Fn = T22) [l^l, is given by [l^ 



where 

A± = Mn±|Mi2|, (21) 

and p is the complex number 

/ = M12. (22) 
From Eq. ((20)) . and in terms of the CP-enantiomers we have 
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|P> = ^(|*+>+|*_ >), (23) 

|i?> = ^(|*+ (24) 

(25) 

where e~*" = p/|p|. 

If on the other hand we assume only T-invariance, then MI2 — M12 [20l |. and in this case the eigenvalues and 
eigenvectors of Eq. are given by [1^1 

A± = ^(Afn + M22) ± \ [(Afii - M22f + AMl^Y'\ (26) 

|*+ > = cos</)|L > +sin0|P >, (27) 
|*_ > = -sin0|L > -hcos^lP >, (28) 

where cot 20 ^ (^^^f > . 
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C. Corrections to the optical activity 

Now we can determine how the inclusion of muhiple states affects the optical activity with respect to the simplest 
HS model. We prepare the state to be initially \L > and assume first only T invariance. Using Eqs. (I^51) - (P5)) . we 
find the wavefunction at any time is given by 

\^{t) (cos2(/)e-*^+* + sin2 0e-*^-')|L > +sin0cos0(e-*^+* - e-'^-*)|i? > . (29) 

Following a sequence of steps similar to those in Sec [III we can evaluate the probabilities , Pr and the optical 
activity in the presence of multi-state transitions. The probabilities to be in state \L > or |i? > at any time t > are 
given by 

P,it) = cos^(At) + ^^"^~f"^% in^(At), (30) 
PR{t) = {-^)sin\At). (3f) 

We find that 

_ \{M,^-M22f+Ml^cos{2M/h) 

where the explicit matrix elements Map are given by Eg. dTTl) and A = ^[(Afii — M22Y + 4M^2]^^^- 
When CPT invariance holds Mn = M22, and then the optical activity is given by 

Q{t) = QmaxCos{2M/n), (33) 

where now, A = |Mi2|. If we denote the oscillation periods TQprp and when CPT or T invariance is imposed, 
then we have the general result that 

'TCPT>'^T, (34) 

so that the oscillation period in the case of the hypothetical CP-enantiomers is longer then that for the P-enantiomers, 
for a given interaction hamiltonian Hi. 

The time-average optical activity for the case of P-enantiomers is 

Omax/t \{Mii - M22Y + Ml^' 

If CPT invariance is assumed, then the time-average of Eq. (p3| is zero: 

f^) (36) 

"max j I 

As mentioned earlier, we write Hi — Sa^ + ec^ + iJ2, where the first two terms act only within the two dimensional 
subspace of the enantiomers \L >, |i? > and H2 induces transitions from these to the other levels |fc >, thus we 
evaluate 

Ma, = m5a,+ < a\5a. + ea.\P > -PP^^ < a^k x k\H2\P > ^ 

k 

If we shut off the perturbation H2 then from Eq. pT|) it is easy to check that Eq. (15^ reduces to the optical activity 
of the isolated two-state system, Eq.®. By the same token, Ea. (|55|) reduces to the Eq.Q in this same limit. 

Note that since M is Hermitian, PL{t) + Pnit) = 1 continues to hold, even allowing for transitions to the other 
states |fc >. If r were not vanishing, then these probabilities would decay exponentially in time. 

IV. FORMAL ANALOGIES TO THE if-MESON SYSTEM 

In the context of symmetry breaking in physics, Wigner pointed out some time ago a strictly formal analogy between 
TT-mesons and chiral molecules [lH. The neutral K'^ meson and its antiparticle X° are related by the combined 
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operations of charge conjugation and parity (CP): jif" >= CP\K^ > [l^]- From this, one defines superpositions 34| 
\Ki -^{\K° > >) and \K2 >= -^{\K° > -\K° >) that are eigenstates of CP: CP\Ki \Ki > and 

CP\K2 —\K2 >■ The chiral molecules are interrelated by the parity operation: \L >— P\R >, \R >~ P\L >, 
and the eigenstates of definite parity are the mixtures |+ '^(\^ ^ 1^" Tf^l^ 

P|+ >= 1+ > and P\— >= — |— >. These algebraic relationships led Wigner to propose a formal analogy between 
neutral kaons and enantiomers, namely, the state-vector associations (o) 

\L> ^ \K" >, (38) 

\R> ^ \K" >, (39) 

1+ > o |i^i >, (40) 

\->^\K2>. (41) 

This analogy can be made more encompassing, by extending these relationships to the dynamic level. The kaons 
are eigenstates of the strong Hst and electromagnetic interactions: {Hst + Hj)\K^ >= mK\K° > and {Hgt + 
Hj)\K'^ >= mif > and are degenerate in mass. The weak interaction H^eak connects K° and if with other 
continuum states which causes the various decay modes and removes their degeneracy. The Schrodinger equation 
Eq.®, describes the time evolution of a neutral kaon system, within the two- level approach in Eas. p^ - (|18p . For 
this, one makes the specific identifications 

Ho = H,t+H^, (42) 

Hi = H^eak, (43) 

(44) 



where now the indices a,^ stand for the states or K'^ and m — rriK — mp^ [l9l - l2li |35| . Due to kaon decays, 
the decay matrix T is nonzero. The eigenvalue problem for (M — iT) has been worked out in full detail [l^. The 
eigenvalues have real and imaginary parts, and these can be expressed as mi, 2 — |7i,2 US]- So if we prepare a state 
which is initially pure at any later time the probability to find a K'^ is '21'] {% — 1) 

P{K°,t) = i[e-''i* + e-'^=*-2e-('^i+^^)*cos(TO2-mi)i], (45) 

= cos^(i!^^^^), (,.2^0), (46) 

which describe decaying meson oscillations. In the limit of zero decay (second line) these would become pure oscil- 
lations and thus formally similar to the chiral oscillations derived above. These considerations raise Wigner's static 
analogy relating state vectors, to a dynamic one, between effective two-level Schrodinger equations, oscillations in the 
transition probability between L and R, Eq. ©, and the strangeness oscillations of the neutral kaon system. 

We must point out however that the consideration given to a two-state model involving CP-enantiomers in Sees 
IIIIBI and IIII CI is strictly a mathematical illustration only and can never be realized experimentally. This is because, 
such molecule-antimolecule transformations would require a huge violation of baryon number conservation [s^l . This 
problem does not arise in the K-meson system because mesons have baryon number zero (mesons are not baryons!) 



V. DISCUSSION 



A number of criticisms of the simplest two-state HS model have been marshalled in the past: namely that it could 
only apply at exceedingly low temperatures, that it neglects the radiation field, and that it does not account for 
collisions. It has been suggested that perhaps the most serious problem arises from the multi-state nature of real 
molecules [l^. It is the latter objection which motivates the work presented here, providing us the incentive to consider 
the multi-state nature in a fairly general way. One aspect (and only one) of this complex problem is the influence 
of transitions to a tower of excited (electronic- vibrational) states of enantiomers induced by an appropriate external 
field or perturbation. In this situation, results from Weisskopf- Wigner perturbation theory demonstrate that we can 
continue to employ an effective two-state description, where the influence of the tower of multiple states is accounted 
for by the matrix elements of the effective hamiltonian or mass matrix acting in the subspace of the two enantiomers. 
The effects that multiple states have on the racemization and optical activity can then be worked out in terms of 
the explicit matrix elements of the speciflc interaction responsible for these transitions. The importance of assuming 
CPT or T-invariance is underscored here. These results hold generally. The main result is that a two-state approach 
remains valid, because the inclusion of the multiple states can be included in an effective Hamiltonian description. 
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An approach such as this may prove useful for interpreting proposed spectroscopic measurements of molecular parity 
violation, such as represented in Fig. 1 of reference [32i] . involving transitions to excited levels. 

In comparison to the general results obtained here, much of the theoretical work has focused on the explicit 
calculations of the parity violating energy difference in chiral molecules; see [s^ for a recent review. These involve 
ab initio computations of the PV interactions employing techniques such as non-relativistic and relativistic (Dirac)- 
Hartree-Fock and multi-configurational self-consistent (MCSCF) levels, as well as density functional theory (DFT) |22| - 
[2p|. There the primary objective is the (numerical) evaluation of effective parity- violating Hamiltonians which requires 
using many-body quantum-mechanical wave functions, to account for the multiple nuclei and electrons ^32i | involved. 
In these investigations, the dominant contribution to the parity- violating energy difference between enantiomers Epy 
is calculated from matrix elements connecting the ground-state singlet with excited triplet states. Thus, it should be 
possible to employ an effective two-level description in these more complex theoretical approaches as well. 

Regarding collisional effects, a two-level approach has been used to describe how racemization depends on the 
interaction of the enantiomers with the environment (3^ . This might suggest that both the multi-state nature of real 
molecules and collisional effects with the surrounding medium might be able to be combined in an overall effective 
two-level description. 

Finally, we have also further developed the formal dynamical analogies between the system of enantiomers and the 
kaon system. The unifying framework is provided by the WW perturbation theory. We note that a formal comparison 
between chiral molecules and neutrinos was recently invoked to derive properties of the oscillations between isolated 
enantiomers in a two-level HS-type approximation |36i] . 
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Appendix A: Reduction from multi-level to a two-state system 



The main steps of the Weisskopf-Wigner (WW) time dependent perturbation theory [17Hl9l | and Appendix A of 
dH, are reviewed and adapted here to a doublet of mass-degenerate enantiomers. We emphasize that the WW method 
is general and provides a way to reduce an a-priori multiple level quantum system to an effective two-level system, 
independent of the actual form of the specific hermitian Hamiltonians involved. 

To derive Eg. p^ from Eqs. p^ - pS]) . introduce the two-component column vectors 

"PW - h(A ^ i,\Tri„^ , (Al) 



6(t) J '^'^ \< k\Hi\R > 

then the first WW approximation consists in truncating the solution of the Eqs. (fT3l) - (fT5)) to second order in Hi. This 
implies they can be written as follows, by using < k\H[\L >= e'"''=* < k\Hi\L >, < k\H[\R >= 6™*=* < k\Hi\R >: 

^ h0(t) + J2 CI-e-™'=*Cfc(i), (A2) 

k 

^dc^ = Cl0(t)e™^*, (A3) 
where Wk = — m, {rriL = mn = m) and h is the submatrix of Hi in the two-state subspace: 



<L\Hi\L> <L\Hi\R> 
< R\Hi\L > < R\Hi\R > 



(A4) 



Solve Ea. (|A3[) for Ck and substitute these solutions back into Ea. (|A2p . The resultant equation for <p{t) can be solved 
in closed form via a Laplace transform [2l| and yields 

1 r°° g(*y+e)t 
'»'"-^L'' %-.. + W(.,H..) ^°- («) 
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where <pQ = 0(0) is the initial condition, and 

W(.)^h-^-^, B.^Cld (A6) 

Up to this point, the solution Eq. ljASP is exact to 0{Hf). If we regard the perturbation Hi as small, then the second 
order contribution to the matrix W should receive its main contribution to the integral from the neighborhood of 
y = 0. The second WW approximation consists in replacing W by its value at y = 0, which leads to the integral 

4>{t) = 7r- dye^2/t y + h-PP^^-i7r^<5K)DJ cj^o, (A7) 

^"■^ J-oo ^ k k 

and which follows from the identity {PP denotes the Cauchy principal part) [31 

1 PP 

lim = =F i7r<5(x). (A8) 

CT-»-o X ±ia X 

Evaluating the integral Eq. (|A7[) yields the general solution 

ct>{t) = e-'^^Vo, (A9) 

where 

Wo = h-PPV^-i7ry<5K)Dfe. (AlO) 

k k 

The time dependence in the interaction representation of the two-level wavefunction in the WW approximation is 
given by 

^±ct>it)^ Worn- (All) 

Returning now to the Schrodinger representation $ = e~'^°*0, 

i^* = {Ho + e-'""'Woe'""'')^, (A12) 
= (i/o+Wo)*, (A13) 
M-ir)*. (A14) 

Here 

M = ml + h-PPV— , (A15) 

and 

r = 27ry <5(w;fc)Dfc, (A16) 

k 

are known as the mass and decay matrices, respectively. Since the two level subsystem is degenerate, Hq = ml and 
so [e-»^«*,Wo] =0. 
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